Synchronization on Effective Networks 



Tao Zhou^'^^ Mij^g zhao^'^'^, and Changsong Zhou=E 
IDepartment of Modern Physics, University of Science and Technology of China, Hefei 230026, PR China 
^Department of Physics, University of Fnbourg, Chemin du Musee 3, CH-1700 Fribourg, Switzerland 
"^Department of Physics, Centre for Nonlinear Studies, 
and The Beijing-Hong Kong-Singapore Joint Centre for Nonlinear and Complex Systems (Hong Kong), 
Hong Kong Baptist University, Kowloon Tong, Hong Kong, China 
"^College of Physics and Technology, Guangxi Normal University, Cuilin 54IOO4, P- R- China 

The study of network synchronization has attracted increasing attention recently. In this paper, 
we strictly define a class of networks, namely effective networks, which are synchronizable and 
orientable networks. We can prove that all the effective networks with the same size have the same 
spectra, and are of the best synchronizability according to the master stability analysis. However, it 
is found that the synchronization time for different effective networks can be quite different. Further 
analysis show that the key ingredient affecting the synchronization time is the maximal depth of an 
effective network: the larger depth results in a longer synchronization time. The secondary factor 
is the number of links. The more links connecting the nodes in the same layer (horizontal links) 
will lead to longer synchronization time, while the increasing number of links connecting nodes in 
neighboring layers (vertical links) will accelerate the synchronization. Our findings provide insights 
into the roles of horizontal and vertical links in synchronizing process, and suggest that the spectral 
analysis is helpful yet insufficient for the understanding of network synchronization. 
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Synchronization is observed in many natural, social, 
physical and biological systems, and has found apphca- 
tions in a variety of fields As a result, the subject 
of synchronization is continuously calling for serious and 
systematic investigation, and has evolved to be an inde- 
pendent field of scientific research. Recently, the dramat- 
ically increasing interests in complex networks have been 
pervading the study of synchronization — understanding 
the synchronizing process of a network of dynamical sys- 
tems is attracting more and more attention (see the re- 
view articles ^ ^ and the references therein). A par- 
ticularly interesting issue is how to enhance the network 
synchronizability Ij]. Under the framework of master 
stability analysis [a if the synchronization region is 
bounded and the nodes are coupled linearly and sym- 
metrically, the network synchronizability can be mea- 
sured by the eigenratio, which is defined as the ratio of 
the largest eigenvalue to the smallest nonzero eigenvalue 
of the coupling matrix. A smaller eigenratio indicates 
a better synchronizability, and vice versa. Some effec- 
tive synchronizability-enhancement methods, aiming at 
reducing the eigenratio, are proposed. These methods 
include the discrete optimization of network structure 
0, p, the regulation of coupling pattern [ p. [Tol . [ill . [H, 
[l3l . [l4l | , the adaptive coupling strategy [la, [la] and the 
structural modification [l3,[ll,[il|. It is worth noting, al- 
though the coupling matrix of coupling pattern regulated 
networks and adaptive networks may become asymmet- 
ric, which correspond to directed networks, the master 
stability analysis is still valid Q. 
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Recently, Nishikawa and Motter [2O, [Hj proposed a 
class of maximally synchronizable networks that have the 
smallest eigenratio, R= 1. They showed that a network 
is maximally synchronizable if (i) it embeds an oriented 
spanning tree, (ii) it does not have any directed loops, 
and (iii) it has normalized input coupling strength of 
each node. Their work [23, [2l[ provides a start point 
to the understanding of the role of directivity and feed- 
back loops in network-based synchronization [221. and 
is also relevant to the synchronization control |23l. [2^ . 
However, some crucial issues, such as the time required 
for full synchronization, lack serious consideration thus 
far. 

In this paper, using the language of graph theory [2^, 
we define a new class of networks, namely effective net- 
works, as the synchronizable and orientable networks, 
and prove that the coupling matrix of any effective net- 
work has eigenvalues = Ai < A2 = ■ ■ • Aat = 1. We 
further propose a preferential selection algorithm (PSA) 
to extract effective network from any given connected 
undirected network, and inve stig ate the synchronization 
of chaotic Logistic systems [2^ on effective networks, 
emphasizing on the relation between the topology and 
synchronization time. Surprisingly, although the spec- 
tra of all effective networks are exactly the same, there 
is a great difference between their synchronization time. 
Based on extensive simulations, it is found that the syn- 
chronization time is mainly determined by the maximal 
depth of an effective network, complementary to a recent 
study about the synchronization of the Kuramoto model 
in an ideal tree [231 ■ The secondary factor affecting the 
synchronization time is the number of links. The more 
the links connecting the nodes in the same layer [hori- 
zontal links), the longer the synchronization time. Our 
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FIG. 1: (Color online) tsync vs. r for diflerent average de- 
grees, {k)—4:, 6, and 10. The dot line indicates the crossover. 
The network size is 1000, and the coupling strength takes 
0.55. Each data point is obtained by averaging over 100 in- 
dependent realizations. 



simulation suggests a logarithmic relation between the 
synchronization time and the number of those links. In 
contrast, the increasing number of links connecting nodes 
in neighboring layers {vertical links) will accelerate the 
synchronization speed. Given the maximal depth, the 
simulation results indicate that when the network con- 
sists of abundant horizontal and vertical links, the ratio 
of the number of horizontal links to the number of verti- 
cal links mainly determines the synchronization time. 

We denote D{V,E) a directed graph, where V is the 
set of node and E the set of directed links. The multiple 
links and self-connections are not allowed. A directed 
network is called synchronizable if for any two differ- 
ent nodes Vi and Vj, 3P{vi Vj), or 3P{vj Vi), or 
3k ^ i,j such that 3P{vk — > Vi) and 3P{vk vj). Here, 
P{vi Vj) denotes an arbitrary directed path from Vi 
to Vj. In language, any pair of nodes are either coupled 
directly or driven by some common nodes. A directed 
network is orientable if it contains no directed cycles. 
Note that, a synchronizable network must be connected, 
while an orientable network is not necessarily connected. 
Wc normalize the input coupling strength of each node: 
hi = 1, hj = —^/k~ if the directed link from Vj to Vi 
exists, otherwise kj = 0. Here, k~ denotes the in-degree 
of node Vi . A node is called source if its in-degree equals 
zero, accordingly, we have a Lemma: Any orientable net- 
work has at least one source. Proof. — Assume there ex- 
ists an orientable network D containing no source. De- 
note wi — > U2 ^ ■ • • — > the longest directed path in 
D. Since k^_^ > 0, there exists at least a link v vi. 
For D contains no directed loops, v ^ V2,--- ,vi, thus 
the directed path v ^ vi ^ ■■■—>■ vi is even longer, 
which conflicts to the above assumption. Therefore, D 
contains at least one source. For effective network (i.e., 
synchronizable and orientable network), we have a The- 
orem: The normalized Laplacian of an effective network, 
L, has eigenvalues = Ai < A2 = ■ • ■ Xn = 1- Proof. — 
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FIG. 2: (Color online) The change of (a) the maximal depth, 
(b) the average depth, (c) the number of vertical links, -/V„, 
and (d) the number of horizontal links, N^, with the parame- 
ter T for average degree (A:)=4, 6, and 10. The network is a BA 
scale-free network with -Af=1000. Each data point is obtained 
by averaging over 100 independent realizations. Symbols are 
of the same meanings as those in Fig. 1. 



There exists an ordered list of V such that every link in 
E is from a former node to a later node. One can get this 
list from an empty list according to the following steps: 
(i) Add all sources of the current network at the top of 
the list in an arbitrary order; (ii) Remove those source 
from the current network and go to step (i) until the net- 
work becomes empty. This process is practicable since 
the network is finite, and after each removal of sources, 
the remain network is still orientable, thus has at least 
one new source (see the Lemma). Denote such an ordered 
list as Q = {vi,V2, ■ ■ ■ , Wat}, clearly. Hi > j, the directed 
link (vi — > Vj) can not exist because in the induced sub- 
graph of nodes {vj, Wj+i, • • • , vn}, Vj is a source without 
any in-link. Shift columns and rows of L corresponding 
to the ordered list Q and denote it by Lq, then, Lq is 
an down-triangle matrix having the same spectrum as L. 
Therefore, all the eigenvalues of Lq, determined by the 
equation det(A/ — Lq) = 0, are its diagonal elements. 
Any node having in-degree larger than corresponds to 
a unit diagonal element, while a source corresponds to 
a zero diagonal element. Clearly, vi is a source, and if 
there exists another source , since there are no directed 
paths end at fi or Ufe, D is not synchronizable, conflicting 
to the condition that D is an effective network. Actually, 
any effective network has one and only one source, so the 
normalized coupling matrix has one zero-eigenvalue and 
— 1 eigenvalues all equal 1 . 

We propose an algorithm, namely preferential selection 
algorithm (PSA), to extract effective networks from any 
given connected undirected network G{V,E). The pro- 
cedure is as follows: (i) Set an empty list Q; (ii) Select a 
node Vi in G with probability proportional to fcf where 
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FIG. 3: (Color online) The synchronization time vs. the max- 
imal depth. The black squares and red circles represent two 
different configurations of the toy model and the former has 
smaller average depth. The inset compares the average depth 
of the two configurations. The network size is A'" = 1001 and 
the number of vertical links is Nv = 1000. Each data point is 
obtained by averaging over 1000 independent realizations. 



ki is the degree of Vi and r is a free parameter. Put this 
node at the top of Q; (iii) Piek a node Vj with probability 
proportional to fcj, satisfying that it does not belong to 
Q and at least one of its neighbors belongs to Q already. 
Push it at the end of Q; (iv) Repeat the step (iii) until 
all nodes in V belong to Q. Since G is connected, the 
termination condition can always be achieved. Denoting 
Q = {vi,V2, ■■■,vn} where A'' is the size of G, then for 
each link {vi ,Vj) G we set its direction as Vi Vj if 
i < or as Vj —> Vi if i > j. All those directed links con- 
stitute a set E', and D{V,E') is an effective network. The 
proof is straightforward and thus is omitted here. Effec- 
tive network has a clear hierarchical structure that the 
nodes closer to the source have higher levels. The PSA 
with positive r tends to put high-degree nodes at higher 
levels. For a given undirected network G, the extracted 
effective networks for different r have far different topolo- 
gies, and we will see below that they exhibit far different 
synchronization behaviors. 

We choose a benchmark individual dynamical system 
for the numerical analysis, namely the Logistic map [1^. 
For an arbitrary node Vi , we denote its state as Xi , which 
evolves as Xi^t+i = f{xi,t 



(BA) model [2^. Given a r, we first extract the effec- 
tive network using the PSA, and then randomly assign 
each node an initial state between and 1, and then all 
the nodes iterate according to the coupling equations for 
2000 time steps, which is long enough to drive the sys- 
tem to the synchronization state. Then, we give each 
node a perturbation, that is, adding a small number ran- 
domly chosen in (-0.005, 0.005). The synchronization 
time, tsync, is defined as the time steps required for the 
system to return the synchronization state. In the numer- 
ical implementation, we monitor the standard deviation, 
CT, of the states of all nodes, and the system is considered 
to be synchronized if cr < 10~^ and will not increase any 
longer. 

Figure 1 reports the synchronization time as a function 
of r. Clearly, the synchronization time decreases remark- 
ably as the increasing of r. Although the effective net- 
works corresponding to different r have exactly the same 
spectrum, their synchronization behaviors are greatly dif- 
ferent and the one with higher-degree nodes in the top- 
per positions synchronizes faster (i.e., a larger r leads 
to a faster synchronizing process). In addition, there is 
a crossover of the curves for different average degrees at 
Tc w (represented by a dot line in Fig. 1). Interestingly, 
when r is larger than Tc, the network with smaller aver- 
age degree will synchronize faster, in contrast, the net- 
work with larger average degree will synchronize faster. 
To understand these phenomena observed in Fig. 1, we 
next study the structure properties of effective networks 
at different t. 

In an effective network, the number of links in the 
shortest path from the source node to a node Vi is called 
the depth of Vi, denoted by d^, and the depth of the 
source node is zero. The maximal depth and the aver- 
age depth are defined 

(d) — X^ili respectively. We classify the links into 
two categories: vertical links are the ones connecting two 
nodes with different depths (obviously, the difference is 
one) while horizontal links connecting two nodes with 
the same depth. Figure 2 reports the change of these 
structural measurements with r. With the increasing of 
T, the number of horizontal links will increase too, while 
the maximal depth, the average depth and the number 
of vertical links will decrease. Especially, the shapes of 
djiax vs. T and {d) vs. r curves are very similar to that 



is the overall coupling strength. Uj is the element of the 
normalized Laplacian L, f{x) = ax{l — a;) is a Logis- 
tic map with a = 4.0 fixed in this paper. According to 
Ref. [2^ , it can be proved that the effective network will 
achieve the complete synchronization state when the cou- 
pling strength satisfies the inequality 1 — < e < 1, 
where A is the Lyapunov exponent of f{x) at the current 
parameter. 

To reveal how the topology affects the efficiency of syn- 
chronization on effective networks, we study the synchro- 
nization time versus the parameter r. The underlying 
networks before PSA are generated by Barabasi- Albert 



s'Y^^ ^ lijf{xjt) where s ^sync vs. r curve, therefore one may guess that either 



c^max or {d) (or both) mainly determines the synchroniza- 
tion time. To judge whether it is true, we propose a toy 
model, where iV -I- 1 nodes are connected by N vertical 
links starting from a single source node, resulting in a tree 
structure. Given the maximal depth dmax (in the present 
simulations, we tune dmax in the range 1 < djiax £ 10), 
we (i) put ^ nodes in each of the dmax ~ 1 layers from 



1 < d < 



1 and the remaining ( 1 



10 



nodes 



in the dmax-th layer; or (ii) put ^ nodes in each of the 
dniax ~ 1 layers from 2 < d < d^ax a-nd the remaining 
(l — '^""o" ^ ) nodes in the first layer. Each d-th layer 
node with 1 < d < dnax connects to (with an in-link) the 
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FIG. 4: (Color online) Effects of (a) the number of horizontal 
links and (b) the number of vertical links on the synchroniza- 
tion time. The network size is = 1001 and each data point 
is obtained by averaging over 1000 independent realizations. 

node in the (d— l)-th layer one by one. Clearly, if d^ax = 1 
or dmax = 10, the two methods lead to statistically the 
same (d), while for 1 < dmax < 10, the former method 
(putting more nodes in the last layer) gives longer (d). 
In the inset of Fig. 3, we report the difference of average 
depths in the case of = 1000. It is clearly observed 
in Fig. 3 that the synchronization time is approximately 
proportional to the maximal depth, while almost not af- 
fected by the average depth. 

Further more, to reveal the roles of the vertical and 
horizonal links, we randomly add some vertical and hori- 
zonal links in the above-mentioned tree topologies. As 
shown in Fig. 4(a), the horizontal links hinder the syn- 
chronization. In the case of minimal number of vertical 
links (i.e., Ny = 1000), the horizontal links, together 
with vertical links, make up some triangles in the form 
{vi Vj,Vi — > Vk,Vj — > Wfe} with Vj and Vk in the same 
layer; and in the cases of abundant vertical links (i.e., 
Ny = 2000 and Ny ~ 3000), they make up two types of 
triangles, respectively with Vj and Vk in the same layer 
and with Vi and Vj in the same layer. These triangles in- 
troduce additional disturbances among individuals, and 
thus slow down the synchronization. The simulations 
suggest a logarithmic relation between <sync and A'^^ at 
large N^. 

The effects of vertical links are twofold. Firstly, more 
vertical links are helpful for the propagation of control- 
ling signals from the upper-layer nodes to the lower-layer 
nodes. At the same time, the additional vertical links will 
increase the number of triangles, and thus hinder the syn- 
chronization. These two competing effects bring nontriv- 



ial role of vertical links. In simulations, we find that when 
there are not so many horizontal links, adding vertical 
links will first hinder synchronization, and then enhance 
synchronization, leading to a maximal synchronization 
time (Fig. 4(b)). When the network possesses abundant 
horizontal links, the more vertical links will accelerate the 
synchronization. Simulation results in Fig. 4(a) and Fig. 
4(b) suggest that when the network consists abundant 
vertical and horizontal links, the ratio Ni-jNy mainly de- 
termines the synchronization time: the larger the ratio, 
the slower the synchronizing process. This gives a qual- 
itative explanation of the crossover observed in Fig. 1. 
When T > 0, the maximal depths of networks with differ- 
ent average degrees are more or less the same, but Nh/Ny 
in the case of (fc) = 10 is much larger than the cases of 
(k) = 4 and (k) = 6. Therefore, when r > 0, a network 
with larger average degree synchronizes even slower than 
the one with smaller average degree. 

In conclusion, we have proposed a preferential selection 
algorithm to extract effective networks from any given 
connected undirected network. Although the spectra of 
all effective networks are the same, their synchronization 
time can be far different. The key ingredient affecting 
the synchronization time is the maximal depth of an ef- 
fective network: the larger the dmax, the longer the tsync- 
The secondary ingredient is the number of links. In the 
cases of abundant links, the smaller number of horizon- 
tal links and larger number of vertical links will lead to 
a faster synchronization. Actually, with fixed c?max, our 
simulation indicates that the ratio N^/Ny mainly deter- 
mines the synchronization time. Under the framework 
of master stability analysis 0, Q, the majority of pre- 
vious studies on network synchronization only focus on 
the spectral analysis of the coupling matrix [2, Q , and 
it has been shown that synchronization time is propor- 
tional to I/A2 for typical networks [lO], but this is not 
the case for the effective networks we propose, which have 
exactly the same spectrum yet exhibit noticeably differ- 
ent synchronization behaviors. Based on the results re- 
ported here, we claim that the spectral analysis is helpful 
but not enough for the understanding of network syn- 
chronization. Our findings not only em pha size again the 
fundamental importance of directivity [20l . [2l| , but also 
provide insights into the roles of horizontal and vertical 
links in the synchronizing process. 
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